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Abstract 

We analyze higher order decay modes of the lightest top squark t\ in the Minimal 
Supersymmetric extension of the Standard Model (MSSM), where the lightest SUSY 
particle (LSP) is assumed to be the neutralino Xi- F° r small t\ masses accessible 
at LEP2 and the Tevatron, we show that the four-body decay mode into the LSP, 
a bottom quark and two massless fermions, t\ — > bxiff, can dominate in a wide 
range of the MSSM parameter space over the loop-induced decay into a charm 
quark and the LSP, i\ — > cxi- This result might affect the experimental searches 
on this particle, since only the later signal has been considered so far. 
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1. Introduction 



Supersymmetric (SUSY) theories, and in particular the Minimal Supersymmetric exten- 
sion of the Standard Model (MSSM) [I], 0, predict the existence of a left- and right 
handed scalar partner, fi and fa, to each Standard Model (SM) fermion /. These cur- 
rent eigenstates mix to form the mass eigenstates f\ and fi- The search for these SUSY 
scalar fermions is one of the main entries of the LEP2 and Tevatron agendas. At the 
Tevatron, the production cross sections of squarks are rather large since they are strongly 
interacting particles and stringent bounds, rriq > 250 GeV |3|, £|, have been set on the 
masses of the scalar partners of the light quarks by the CDF and DO collaborations. At 
LEP2, bounds close to the kinematical limits, m\ > 80 GeV [|, have been set on the 
masses of the charged scalar leptons, while the experimental bound on the mass of the 
sneutrinos is still rather low, ma > 45 GeV H . 

The situation of the top squarks is rather special. Indeed, the two current stop eigen- 
states and tn could strongly mix |J due to the large m t value which enters in the 
non-diagonal element of the mass matrix. This leads to a mass eigenstate t\ possibly 
much lighter than the other squarks, and even lighter than the top quark itself. If the 
stop t\ is lighter than the top quark and the chargino [and also lighter than the scalar lep- 
tons], the two-body decay modes JFj into a top quark and the lightest neutralino [which, 
in the MSSM with conserved R-parity 0, is expected to be the lightest SUSY particle 
(LSP)] and into a bottom quark and the lightest chargino, are kinematically forbidden at 
the tree-level. The main t\ decay channel is then expected to be the loop-induced and 
flavor-changing decay into a charm quark and the lightest neutralino [[| 

ti -> c X ? (1) 

At the Tevatron, a light scalar top squark can be produced either directly in pairs 
through gluon-gluon fusion and quark-antiquark annihilation, gg/qq — > tit* [10|, or in 
top quark decays, t — > t\\i [PL if kinematically allowed. With the assumption that the 
branching ratio for the decay t% — > CX\ is 100%, a contour in the m^-m x o plane has been 
excluded by the CDF collaboration in a preliminary analysis [12 ]. For instance, for a 
neutralino mass of m x o ~ 40 GeV, the maximum excluded t\ mass is m~ tl ~ 120 GeV; 
for smaller or larger m x o values the bounds on are lower. At LEP2, the lightest top 
squark is produced in pairs through s-channel photon and Z-boson exchange diagrams, 
e + e~ — > 7, Z — > t\t\ |Tj|. Again, assuming a branching ratio of 100% for the decay 



t\ — > cXi, the LEP collaborations have set a lower boundQ of > 83 GeV on the 
lightest stop mass, with the additional assumption that the amount of missing energy is 
larger than 15 GeV 0. 

1 Recently the OPAL collaboration [Q has set a stronger bound of > 87.2 GeV with a mass 
splitting between the stop and the LSP larger than 10 GeV. 
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All these searches rely on the fact that the decay t\ — > ex? is largely dominant^. 
However, there is another decay mode which is possible in the MSSM, even if the lightest 
top squark is the next-to-lightest SUSY particle: the four-body decay into a bottom 
quark, the LSP and two massless fermions 

h - b X °Jf' (2) 

This decay mode is mediated by virtual top quark, chargino, sbottom, slepton and 
first/second generation squark exchange, Fig. 1, and is of the same order of perturba- 
tion theory as the loop induced decay t\ — ► c%°, i.e. 0(a 3 ). In principle, it can therefore 
compete with the latter decay channel. Several estimates of the order of magnitude of the 



decay rate of the process eq. (2) have been made in the literature || [LJJ . These estimates 
were based on the assumption of the dominance of one of the contributing diagrams [the 
last diagram of Fig. lc in Ref . for instance] , and the exchanged particles are assumed 
to be much heavier than the decaying stop squark [therefore working in the point-like 
limit to evaluate the amplitudes]. In this case, the output [as one might expect since the 
virtual particles were too heavy from the beginning] was that the decay rate eq. (2) is in 
general much smaller than the decay rate of the loop induced decay into a charm quark 
and a neutralino. 

The purpose of this paper is to revisit the four-body decay eq. (2) in the light of the 
recent experimental limits on the masses of the SUSY particles. We perform a complete 
calculation of the decay process, taking into account all Feynman diagrams and interfer- 
ence terms. We show that, if the exchanged particles do not have a too large virtuality [i.e. 
that they are not much heavier than the decaying top squark] this four-body decay can 
in fact dominate over the loop induced decay channel eq. (1) in large areas of the MSSM 
parameter space. This result will therefore affect the present experimental lower bounds 
on the ti mass, since as discussed previously, this state has been searched for under the 
assumption that the decay mode t\ —> c\\ is the main decay channel. A fortran code 



calculating the partial widths and branching ratios for this four-body decay mode |T7] 



is made available, and the lengthy formulae for the four-body decay width will be given 



elsewhere [17 



The rest of the paper is organized as follows. In section 2, we introduce our notation 
and discuss the two-body decay modes of the top squarks, and in particular the decay of 
the lightest stop squark into charm and neutralino, paying special attention to the cases 
where the decay rate can suppressed. In section 3, we analyze the four-body decay mode 
eq. (2), and make a detailed numerical comparison with the previous decay channel. Some 
conclusions are then given in section 4. 



2 Some three body decays of the top squark have also been considered [[L5). At LEP2 for instance, the 
possibility of a light sneutrino, leading to the kinematically accessible decay mode t\ — ► blu, has been 
analyzed; in this case, a bound > 85 GeV Q has been set on the lightest stop mass. 
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2. The Two-Body Decays 



In this section, we will first summarize the properties of top squarks: masses and mixing, 
and then discuss their tree-level two-body decays into neutralinos and top quarks, and 
charginos and bottom quarks as well as the loop induced decay of the lightest top squark 
into a charm quark and the lightest neutralino. 

2.1 Squark masses and mixing 

As discussed previously, the left-handed and right-handed squarks of the third generation 
fi and fn [the current eigenstates] can strongly mix to form the mass eigenstates f\ and 
/ 2 ; the mass matrices which determine the mixing is given by 



M) = 



m LL m f A f 

m f A f m 2 RR 



(3) 



with, in terms of the soft SUSY-breaking scalar masses m? and m? , the trilinear cou- 
pling Af, the higgsino mass parameter /i and tan/5 = vu/vd, the ratio of the vacuum 
expectation values of the two-Higgs doublet fields 

m 2 LL = m 2 f + m 2 j L + (if - e f s 2 w ) cos 2/3 Mf 
m RR = m / + m / fl + e f s w cos 2/3 M 2 Z 
A, = A f -^npy 2I f (4) 

with ef and 7| the electric charge and weak isospin of the sfermion / and s^y = 1 — = 
sin 2 6 W . The mass matrices are diagonalized by 2 x 2 rotation matrices of angle 9f 





VJ — \ 6f 6f , c 8f =cos9 f and s$ =sm9f 
V - s 0f c °f J 

The mixing angle 9f and the squark eigenstate masses are then given by 
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(5) 



(6) 



(7) 



In our analysis, we will take into into account the mixing not only in the stop sector 
where it is very important because of the large value of m t , but also in the sbottom and 



stau sectors, where it might be significant for large values of the parameters \i and tan (3. 
Furthermore, we will concentrate on two scenarii to illustrate our numerical results: 

(i) "Unconstrained" MSSM: we will assume for simplicity a common soft SUSY- 
breaking scalar mass for the three generations of squarks and sleptons and for isospin up 
and down type particles: m- th = m^ R = m~ bR = rriq and nif L = nif R = mj [m~ bh =mi L and 
my L =irtf L by virtue of SU(2) invariance]. The splitting between different particles, and 
in particular between the two top squarks, will be then only due to the D-terms and to 
the off-diagonal entries in the sfermion mass matrices. [Note that in this case, in most 
of the parameter space, the stop mixing angle is either close to n/2 (no mixing) or to 
±7r/4 (maximal mixing) for respectively, small and large values of the off-diagonal entry 
m t A t of the matrix M~.\ Furthermore, we will assume that the mixing between different 
generations is absent at the tree level [otherwise the decay mode t\ — > tXi will occur 
already at this level]. 

(ii) Constrained MSSM: or the minimal Supergravity model (mSUGRA) where the 
soft SUSY breaking scalar masses, gaugino masses and trilinear couplings are universal 
at the GUT scale; the left- and right-handed sfermion masses are then given in terms of 
the gaugino mass parameter mi/2, the universal scalar mass mo, the universal trilinear 
coupling A and tan /3. For the soft SUSY breaking scalar masses, the parameter /i and 
the trilinear couplings at the low energy scale, we will use the approximate formulae for 
the one-loop RGEs given in Ref. fl8| . In mSUGRA [in the small tan/? regime], due to 
the running of the (large) top Yukawa coupling, the two stop squarks can be much lighter 
than the other squarks, and in contrast with the first two generations one has generically 
a sizeable splitting between m~ L and m~^ at the electroweak scale. Thus, even without 
large mixing t\ can be much lighter than the other squarks in this scenario. 

2.2 Two— body decays 

If the top squarks are heavy enough, their main decay modes will be into top quarks and 
neutralinos, U — > tXj [j=l— 4] , and bottom quarks and charginos, U — > bxf [j=l— 2] . The 
partial decay widths are given at the tree-level by 

m 2 x o) - Aa\ j b\ j mim x oe Xj 



m 2 x +) - 4l t ij k t ij mbm x + 

(8) 

where X(x, y,z) = x 2 + y 2 + z 2 — 2 (xy + xz + yz) is the usual two-body phase space 
function and e Xj is the sign of the eigenvalue of the neutralino Xj- The couplings a 1 ^ ^ 
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and b l L j R j for the neutral decay are given by 

m t N j4 



with 



°y 
ft* 



\/2M w sin (3 

m t N j4 
V2M W sin (3 




(9) 



v/2 
-V2 



2 Ar , /l 2 2 

3c w J 



^2 



^1 



(10) 



while the couplings l^j^j and fc^,- for the charged decay mode are given by 



ut 

fi-y 

It 

l 2j 




11' 



In these equations, A" and [//V are the diagonalizing matrices for the neutralino and 
chargino states [19| with 



N' = c w N n + s w N : 



J'2 



-s w N n + c w N i2 . 



(12) 



If these modes are kinematically not accessible, the lightest top squark can decay into 
a charm quark and the lightest neutralino. This decay mode is mediated by one-loop 
diagrams: vertex diagrams as well as squark and quark self-energy diagrams; bottom 
squarks, charginos, charged W and Higgs bosons are running in the loops. The flavor 
transition b — > c occurs through the charged currents. Adding the various contributions, 
a divergence is left out which must be subtracted by adding a counterterm to the scalar 
self-mass diagrams. In Ref. ||, it has been chosen to work in the minimal Supergrav- 
ity framework where the squark masses are unified at the GUT scale; the divergence is 
then subtracted using a soft-counterterm at Agut, generating a large residual logarithm 
log(A GUT /M^/) in the amplitude when the renormalization is performed. This logarithm 
gives the leading contribution to the ti — > cxI amplitude^. 

Neglecting the non-leading (constant) terms as well as the charm-quark mass, the 
partial width of the decay t\ — > c\i is given by 



a 



rn 



LI 



(13) 



3 In fact, this contribution is due to the stop-charm mixing induced at the one loop level, and which 
can be obtained by solving the renormalization group equations in mSUGRA [pf. 
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where fn is given by 



f fc\ 2 ( AT AT 2 2 \ SwNu + C W N 12 

fti = V 2 -{c w N n + s w N 12 ) + - - - s w ) 

.3 \2 3 J CwSw 

and e denotes the amount of Il,r-Cl mixing and is given by ||: 

A L Cg t - A R Sg t 



(14) 



2 2 

m~ r — mf 



(15) 



with 



Al = — V log ( ] ' tb ' co,, : b [m\ +ml +m 2 H + A 2 b 
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^s 2 w 




2M^ cos 2 /? 



47rsfr V il% / 2M^ cos 2 /3 

Assuming proper Electroweak symmetry breaking, the Higgs scalar mass can be 
written in terms of /x, tan /3 and the pseudoscalar Higgs boson mass Ma as 

m 2 Hi = M\ sin 2 (3 - cos 2(3M^ - p 2 (17) 

Note also that and A# are suppressed by the CKM matrix element V c b ~ 0.05 and 
the (running) b quark mass squared m 2 ~ (3 GeV) 2 , but very strongly enhanced by the 
term log (A GUT /M^) which is close to ~ 65 for A GUT ^ 2 ■ 10 16 GeV. 

Let us now discuss the various scenarii in which the decay rate eq. (13) is small. 

(z) First, and as discussed previously, the large logarithm log (Aq UT /M^) ~ 65 appears 
only because the choice of the renormalization condition is made at the GUT scale. This 
might be justified in the framework of the mSUGRA model, but in a general MSSM where 
the squark masses are not unified at some very high scale such as A GU t, one could chose a 
low energy counterterm; in this case no large logarithm would appear. In fact, one could 
have simply made the renormalization in the MS (or DR) scheme, where the divergence 
is simply subtracted, and we would have been left only with the (very small) subleading 
terms. 

(ii) If the lightest top squark is a pure right-handed state [a situation which is in 



fact favored by the stringent constrains [2D| from high-precision electroweak data, and in 
particular from the p parameter; see Ref. pi] for instance], there is no mixing in the stop 
sector, and the e term in eq. (15) involves only the Ar component. For moderate values 
of the trilinear coupling Ab, this component can be made small enough to suppress the 
decay rate T(ti — ► cxD- In addition, the charm squark mass can be made different from 
the lightest stop mass, and taken to be very large; there will be then a further strong 
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suppression from the denominator of eq. (15) [this situation occurs in fact also in the 
mSUGRA model, since because of the running of the top Yukawa coupling, Ir can be 
much lighter than cl, especially for large values of the parameter m ; see Ref. [IBfll- 

(Hi) Even in the case of mixing, for a given choice of the MSSM parameters, large 
cancellations can occur between the various terms in the numerator of eq. (15). Indeed, 
for some values of the soft SUSY-breaking scalar masses and trilinear coupling Aj,, the 
two terms Al,r weighted by the sine and cosine of the mixing angle might cancel each 
other; this would happen for a value of 9 t such that tan# t ~ A^/A/j. In addition, the 
coefficient fi in eq. (14), which summarizes the gaugino-higgsino texture of the lightest 
neutralino might also be very small. In fact, as can be seen, the parameter involves only 
the gaugino components Nu and iV 12 ; if the neutralino Xi is higgsino-like, these two 
components are very small, leading to a very small value for fi. 

Thus, there are many situations in which the decay rate T(ti — > cxi) might be very 
small, opening the possibility for the four-body decay mode, to which we turn our atten- 
tion now, to dominate. 



3. The Four— Body decay mode 
3.1 Analytical Results 

The four-body decay mode t\ — > bxlff, which occurs for stop masses larger than m& + 
m x i, proceeds through several diagrams as shown in Fig. 1. There are first the iy-boson 
exchange diagrams with a virtual top quark, bottom squark or the two charginos states 
[Fig. la]. A similar set of diagrams is obtained by replacing the VT-boson by the charged 
Higgs boson H + [Fig. lb]. A third type of diagrams consist of up and down type slepton 
or first /second generation squark exchanges [Fig. lc]. 

We have calculated the amplitude squared of the decay mode, taking into account 
all these diagrams and interferences. The complete expressions are too lengthy and will 
be given elsewhere JL7|| . We have taken into account the 6-quark mass [which might be 
important for nearly degenerate stop and LSP masses] and the full mixing in the third 
generation sector. We have then integrated the amplitude squared over the four-body 
phase space using the Monte-Carlo routine Rambo |22[ , to obtain the partial decay width. 
Let us summarize the main features of the result. 

The charged Higgs boson exchange diagrams Fig. lb: 

These diagrams do not give rise to large contributions for two reasons. First, because 
of the relation in the MSSM between the charged and pseudoscalar Higgs boson masses, 
M\± = M\ + M^, and the experimental bound Ma > 90 GeV, the charged Higgs boson 
with mass a M H ± > 120 GeV has a much larger virtuality than the W boson contribution; 
since the other exchanged particles are the same, the H contribution is much smaller 
than the W contribution. In addition, the contributions are suppressed by the very tiny 
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Yukawa couplings of the H bosons to leptons and light quarks, except in the case of the 
H + vr coupling which can be enhanced for large tan/3 values; however, in this case the 
decay width T(ti — > cxI) [which grows as 1/ cos 2 /?] is also strongly enhanced. Therefore, 
these contributions can be safely neglected in most of the parameter space. 

The squark exchange diagrams in Fig. la and Fig. lc: 

In general, the diagrams of Fig. lc give very small contributions for t\ masses of the 
order of 100 GeV, since the first and second generation squarks are expected to be much 
heavier, trig > 250 GeV, and their virtuality is therefore too large. For much larger t\ 
masses, the two-body decay channel t\ — > bxt [or one of the three body decay modes] 
will be in general open and will largely dominate the rate. The sbottom contribution 
in Fig. la is also very small, if the mixing in the sbottom sector is neglected and the 
Tevatron bound m~ h > 250 GeV applies [this bound is valid only if bottom squarks are 
approximately degenerate with first /second generation squarks; in the general case the 
experimental bound p| is lower] . For large values of tan (3 which would lead to a strong 
mixing in the sbottom sector with a rather light b\, the decay width t\ — > c\\ becomes 
very large as discussed previously, leaving little chance to the four-body decay mode to 
occur. 

Top quark exchange diagrams Fig. la: 

The contribution of the diagram with an exchanged top quark is only important if the stop 
mass is of the order of m t + m x o and therefore > 0(250 GeV), with the proviso that: 
(i) the lightest chargino must be heavier than t\ for the two-body decay mode t\ — > bxt to 
be kinematically forbidden, and (ii) the t± mass must be smaller than m x o + M\y to forbid 
the three body decay t\ — > bx\W. In models with gaugino mass unification, M 2 ~ Mi/2, 
the Xi an d Xi masses are related in such a way that the above conditions [with a not 
too virtual top quark] are fulfilled only in a marginal area of the parameter space. This 
contribution can be thus also neglected. 

Slepton exchange diagrams Fig. lc: 

In contrast to squarks, slepton [and especially sneutrino] exchange diagrams might give 
substantial contributions, since sleptons masses of (9(100 GeV) are still experimentally 
allowed. In fact, when the difference between the lightest stop, the lightest chargino and 
the slepton masses is not large, the diagrams Fig. lc will give the dominant contribution 
to the four-body decay mode, with a rate possibly much larger than the rate for the loop 
induced decay t\ — » c%°, for small enough values of tan/?. 

Chargino exchange diagram Fig. la: 

The most significant contributions to the four-body decay mode will come in general from 
this diagram, when the virtuality of the chargino is not too large. In particular, for an 
exchanged xf with a mass not much larger than the experimental lower bound m + > 95 
GeV [but with m x + < m x o + Mw to forbid the three-body decay mode], the decay width 
can be substantial even for top squark masses of the order of 80 GeV. 



9 



Thus, a good approximation [especially for a light top squark ~ (9(100 GeV)] is 
to take into account only the lightest chargino and slepton exchange contributions. In 
terms of the momenta of the various particles involved in the process [all momenta are 
outgoing, except for the momentum of the decaying stop which is ingoing] and defining 
the propagators as Dx = Px ~ ^xi the amplitude squared for the charged slepton and 
sneutrino contributions is given by [note that m x o is the physical LSP mass and e Xl is the 
sign of the eigenvalue of the neutralino Xi] 



\M 



4e 6 



„6 n2 7")2 



rnl+a?(U 11 )\l{ 1 ) 2 (p l .ptf)( Pb .p v ) 



(18) 
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-e Xl m x om x + 



(Pb-Pv)(pi-P x +) - (Pb-Pl)(Pu-P x +) + (Pb-P x +)(PvPl) 



where we have taken into account only the exchange of the lightest chargino, and neglected 
the mixing in the charged slepton sector [which is in general not too large for small values of 
tan (3\. In terms of the elements of the matrices diagonalizing the chargino and neutralino 
mass matrices, the couplings a/ and l\ ± read 



l\i = ~ce t Vu + s 6t Y t V 12 , Y t 



m t 



V2M W sin f3 
N 12 



a f = -V2e f N u s w - V2(l! o - e f s 2 w ) 



cw 



(19) 



Taking into account only the contribution of the lightest chargino, the amplitude squared 
for the chargino exchange diagram is given by: 
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(20) 



with N c the color factor, and the chargino-neutralino-W couplings O u ' [the coupling r n 
was given above] 
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(21) 
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If the W boson in the virtual chargino diagram decays leptonically, one has also to take into 
account the interference between the chargino and slepton exchange diagrams. Assuming 
again no mixing in the charged slepton sector [which means that only left-handed sleptons 
contribute], the interference term is given by: 



Re[A xt Ai 



s e w V2Dl + D w D l 



[20 u 2(p b .p xt )(p x +. Pl )(p v .p x o) - (p b .pi)p 2 x +(pv.p x o 



X\> 



x, m xi 



m x\ ll [(Pb-Pl)(P x +-Pu) - (Pb-Pu)(P x +-Pl) - (Pb-P x +)(PvPl) 



4e b (Z f 11 ) 2 a^ 11 
' s 6 w V2Dl + D w D: 



xr 



(22) 



+e xi m x o0 11 [(p b .p v )(p x +.pi) + ( Pb -P x +)(PvPi) ~ (Pb-Pi)(p x +-Pu)} 



Xi 



Xi 



xi 



Note that if the exchanged W bosons and sleptons are real [i.e. for t\ masses larger than 
Mw + m x o + m b and/or mi + m b respectively], the three body decay channels i\ — > bW + X\ 
and t\ — > bll' open up. This situation can be handled by including the total widths of the 
W boson and sleptons in their respective propagators. In the numerical analysis though, 
we will concentrate only on the kinematical regions where the W boson and sleptons are 
off mass-shell. 



3.2 Numerical Results 

For the numerical results, where we include the contributions of all diagrams [and not 
only the dominant chargino and slepton exchange diagrams discussed above] we first show 
in Figs. 2 and 3, the branching ratios for the four-body decay mode BR(t! — > bxiff) 
in the unconstrained MSSM scenario, with a common squark mass rriq. For the gaugino 
sector, we have chosen tan/3 = 2.5, a gaugino mass parameter M 2 = 120 GeV [Figs. 2a 
and 3a] and 200 GeV [Figs. 2b and 3b] and three values of the parameter \x = 300,450 
and 700 GeV. This leads to the lightest chargino and neutralino masses shown in Table 1; 
these values allow for charginos xf with a not too large virtuality, but still heavy enough 
to comply with the available experimental bounds [even for the choice M 2 = 120 GeV]. 
The trilinear coupling A t is varied to fix m~ tx to a constant value, while the coupling A b 
is fixed to A b = —100 GeV. For the additional parameters which enter the t\ — > c\\ 
amplitude, we will take [for the entire numerical analysis, i.e. also for the other figures]: 
Ma = 500 GeV [except in the mSUGRA scenario discussed later, where M4 is given by 
the RGE's], V cb = 0.05, m T b nn - = 3 GeV; the cut-off A will be taken to be the GUT scale 
A GU t = 2 • 10 16 GeV. 

In Fig. 2, the branching ratio BR(ti — > bxiff) is shown as a function of the soft 
SUSY breaking scalar mass trig for two values of the lighter top squark mass = 80 
GeV [2a] and 150 GeV [2b]. The common slepton mass is taken to be mj ~ 200 GeV 
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M 2 


fi 


m v o 


m xt 


80 


120 


300 


50.2 


96 






450 


53.6 


106 






700 


55.4 


112 


150 


200 


300 


87.8 


163 






450 


91.8 


181 






700 


93.8 


190 



Table 1: Xi an d Xi masses for tan/3 = 2.5 and the choices of Mi and fi parameters of 
Figs. 2-3. All parameters and masses are in GeV. 

[so that the contribution of the sleptons to the decay will be extremely small especially 
for the small values]. One sees that even for a rather light stop, m- tx = 80 GeV 
[Fig. 2a], the branching ratio can be largely dominating. For fi = 300 GeV, it is already 
the case for values of rriq close to ~ 200 GeV [which is needed to keep the masses of the 
first and second generation squarks larger than the present experimental lower bound]; 
for larger rriq values, the branching ratio becomes very close to one. For the value fi = 450 
GeV, the branching ratio drops to the level of 10 to 20%; this is mainly due to the 
positive interference generated by the soft scalar mass squared m 2 Hi ~ —fi 2 [eq. (17)] 
which enhances the decay rate T(t\ — > cxi), but also to the fact that the exchanged 
chargino [which gives the largest contribution to the four-body decay channel] is heavier 
than in the previous case, so that its large virtuality suppresses the four-body decay 
mode. For the value fi = 700 GeV, the interference can be either positive or negative, 
and in the margin rriq ~ 300-400 GeV, the decay ii — > &%?//' dominates over the loop 
induced decay i\ — > cxi in spite of the even larger virtuality of the chargino. 

For a heavier t\ state, = 150 GeV [Fig. 2b], and for the value Mi = 200 GeV 
[leading to heavier charginos and neutralinos] and the same fi values as above, the branch- 
ing ratio BR(ti — > bxiff) is larger than 95% for almost all values of the parameter rriq. 
This is first due to the fact that the phase space is more favorable in this case [i.e. the 
virtuality of the chargino is relatively smaller], and also to the fact that for this particular 
choice of the lightest stop mass, the parameter e in eq. (15) which governs the magnitude 
of the decay rate T(t\ — > cxi) is suppressed [i.e. the mixing angle is such that, there is a 
partial cancellation between the two terms in the numerator of eq. (15)]. 

Figs. 3 show the branching ratio BR(tx — > bxiff) as a function of the lightest stop 
mass, for a fixed value of the common squark and slepton mass rriq = mj = 400 GeV, 
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and for the same parameters tan/3, /i and M 2 as in Fig. 2. As can be seen, the branching 
ratio is very small in the lower stop mass range where the virtuality of the exchanged 
chargino is rather large, and increases with increasing m~ tx to reach values close to unity 
near the mi ~ mj, + m + threshold where the two-body decay mode t\ — > bxt opens 
up [of course, beyond this threshold, the loop induced decay becomes irrelevant and we 
stopped the curves at these values]. Note that even for the small values ~ 80 GeV 
[Fig. 3a], the branching ratio for the four-body decay mode ti — > bxiff can reach the 
level of 90%. Since top squarks with these mass values have been experimentally ruled 
out by LEP2 [and possibly Tevatron] searches under the assumption that they decay most 
of time into charm quarks and the lightest neutralinos, the searches at LEP2 have to be 
reconsidered in the light of the possible dominance of the four-body decay mode. 

In the previous figures, sleptons were too heavy to contribute substantially to the decay 
rate t\ — > bxiff since we assumed the common slepton mass to be m z ~ = rriq > 200 GeV. 
In Fig. 4, we relax the assumption m,f = rriq and show the branching ratio BR(ti — > bxiff) 
as a function of the sneutrino mass mo [the masses of the other slepton are then fixed and 
are of the same order] for m~ tx = 80 GeV, /i = 300 GeV [Fig. 4a] and 700 GeV [Fig. 4b] 
and three values of the common soft scalar quark mass m^ = 300, 500 and 800 GeV. As 
can be seen, the contribution of sleptons can substantially enhance the four-body decay 
branching for relatively small masses [corresponding to m„ < 120 GeV in this case]. For 
larger sneutrino masses, the sleptons become too virtual and we are left only with the 
contribution of the lightest chargino discussed previously [and which is constant in this 
case] . 

Finally, let us turn to the case of the mSUGRA scenario. The branching ratio of the 
decay t\ — > bxiff for a light top squark m~ tx ~ 70-130 GeV is shown in Fig. 5, as a 
function of m^ for tan (3 = 2.5, /j, > [Fig. 5a] and // < [Fig. 5b] and several choices of 
the value of the gaugino mass parameter m±/2 ~ O.8M2 [again, the choice of mi/2 leads to 
chargino masses not too much larger than the allowed experimental bounds, m + < 140 
GeV; in fact in this scenario, /i is always large and the neutralinos and charginos and 
almost bino and wino like, with masses m x o ~ M 2 /2 and m x + ~ M 2 ]. Again, one sees 

that BR(t! — > bxiff) can be very large, exceeding in some cases the 50% level, even for 
values mi x ~ 80 GeV, which are experimentally excluded by the negative search of the 
h c Xi signature if this decay channel dominates. 

4. Conclusions 

We have analyzed the four-body decay mode of the lightest top squark into the lightest 
neutralino, a bottom quark and two massless fermions, t\ — > bxiff, in the framework of 
the minimal supersymmetric extention of the Standard Model, where the neutralino Xi is 
expected to be the lightest SUSY particle. Although we have evaluated the partial decay 
width taking into account all the contributing diagrams [and their interferences] , we have 
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singled out those which give the dominant contributions. 



For small t\ masses accessible at LEP2 and the Tevatron, we have shown that this 
four-body decay mode can dominate over the loop-induced decay into a charm quark 
and the LSP, t\ — > CX\, if charginos and sleptons have masses not too much larger than 
their present experimental bounds. This holds in the case of both the "unconstrained" 
and constrained (mSUGRA) MSSM. 

This result will affect the experimental searches of the lightest top squark at LEP2 and 
at the Tevatron, since only the charm plus lightest neutralino signal has been considered 
so far in these experiments. However, the topology of the four-body decay is similar to 
the ones of the three body decay mode t\ — > bll' [for final state leptons] which has been 
searched for at LEP2 || and of the two-body decay mode t\ — > bxf which has been 



looked for at the Tevatron [T3]. The extension of the experimental searches to the decay 



mode ti — > bxiff should be thus straightforward. 
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Figure 2: The branching ratio BR^x — > bxiff) as a function of the common squark mass 
rriq for a scalar top mass m~ tl = 80 (a) and 150 GeV (b). 
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Figure 3: The branching ratio BR(tx — > bxiff) as a function of the top squark mass 
for a scalar mass rriq = 400 GeV and a gaugino mass M 2 of 120 (a) and 200 GeV (b). 
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Figure 4: The branching ratio BR(t! — > bxiff) as a function of the sneutrino mass mi 
for a scalar top mass = 80 and two values of ji — 700 GeV (a) and 300 GeV (b). 
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Figure 5: The branching ratio BR(tx — > bxlff) as a function of stop mass in the 
mSUGRA scenario for ji < (a) and // > (b). 
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